Euler Walk on a Cay ley Tree. 

A.E. Patrick^ 



Abstract. We describe two possible regimes (dynamic phases) of the Euler walk 
on a Cayley tree: a condensed phase and a low-density phase. In the condensed 
phase the area of visited sites grows as a compact domain. In the low-density phase 
the proportion of visited sites decreases rapidly from one generation of the tree to 
the next. We describe in detail returns of the walker to the root and growth of 
the domain of visited sites in the condensed phase. We also investigate the critical 
behaviour of the model on the line separating the two regimes. 

KEY WORDS: Branching processes; critical exponents; martingales; random 
walks. 



1 Introduction. 

Consider a Cayley tree with arrows attached to every site. Initially the arrows point 
at one of the adjacent sites randomly and independently of each other, see Fig. 1. 
An Eulerian walker moves over the Cayley tree according to the following rules. At 
time instants I = 0, 1, 2, . . . the walker jumps from its current location x(l) (at one 
of the sites of the tree) to the adjacent site in the direction of the arrow at x(l). 
At the time of jump the arrow at x(l) is rotated clockwise, till it points to another 
adjacent site. 

We assume that the generations of the tree are numbered from bottom to top. 
Zero generation of the tree contains only the root. If k > I, then the generation 
number k is above the generation number I on a picture of the tree, and we say that 
the k th generation is higher than the generation number /. 

Eulerian walkers were introduced by Priezzhev et al [10] (see [9] for further 
investigations) as a model of "self-organized criticality" . In their version of the 
model the walker rotates arrows on arrival at a site of the graph. The two versions 
are largely equivalent, but in the version used in this paper it is easier to see on the 
picture of a graph where the walker actually goes over the next few steps. 

An attractive feature of the model on a finite graph with reflecting boundary 
(the boundary arrows point inside the graph) is that eventually the walker settles 
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Figure 1: A Cayley tree with internal arrows arranged at random and independently 
of each other. Boundary arrows always point inside the tree to prevent the walker 
from falling off the tree. At site a the walker jumps in the direction of the arrow 
to site b, and the arrow at a is rotated clockwise to point at site c. After that the 
walker jumps back to site a and then jumps to site c. At the time of the last jump 
the arrow at a is again rotated clockwise to point at site d. 



into an Euler circuit, where it passes every edge of the graph twice (once in every 
direction). There is exactly one "clockwise" Euler circuit for any tree, see Fig. 2. 
As a result of the walker's activity, initially chaotically oriented arrows arrange into 
an organized configuration directing the walker around that circuit. 

If the Euler circuit on a tree is a critical state is not that clear. It is not difficult 
to calculate correlation functions for orientations of the arrows at two sites of a 
Cayley tree, assuming the uniform distribution of the current location of the walker. 
For instance, let a\ and a 2 be two arrows at sites in the generations k and fc + mof 
a finite tree containing n generations in total. Then 

Pr[ai =\, a 2 =\] - Pr[oi =\] Pr[a 2 =\] -> -2- 2k ~ m , as n - oo. 

Hence, we have an exponential decay of correlations with the distance between the 
arrows measured in generations of the tree. The same asymptotic behaviour we 
obtain for all other correlation functions, although some of those are positive. 

The above decay of correlations is in contrast to the behaviour found by Dhar 
and Majumdar for the self-organized state of a sand pile on a Cayley tree, see [I]. 
Dhar and Majumdar found that the correlation functions decay as 4 _m , where m is 
the distance between the two sites of the tree. They concluded that the correlations 
are short-ranged, because even after multiplication by the branching factor 2 m one 
still has an exponential decay to 0. Nevertheless they classified the self-organized 
state of the sand pile on a Cayley tree as critical, presumably, because of power-law 
tails in the distributions of avalanche- sizes and related quantities. 
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Figure 2: The Euler circuit on a Cayley tree. The orientation of arrows corresponds 
to the current position of the walker at the root of the tree. 

In our case the correlations decay as 2 _m , and do not vanish after multiplication 
by the branching factor. Therefore one certainly can not rule out the criticality of 
the Euler circuit on a tree on the basis of exponential decay of correlation functions. 
Nevertheless one feels that the self-organized state in this case is closer to the min- 
imally stable state of the ID sand pile, described in the paper j2], than to a truly 
critical state. 

Our main goal in this paper is a description of the formation of an organized 
structure on an infinite tree. We will show that, unlike what one sees on finite graphs 
with reflecting boundary, on an infinite tree a (substantial density of) organized 
structure is not always formed. Of course, if an organized structure is not formed 
on an infinite tree, it is highly sensitive to the boundary conditions and appears on 
a finite tree only as a result of numerous bounces of the walker against the reflecting 
boundary. 

To set the scenery for the study of Euler walk let us describe two possible regimes 
of evolution: a condensed phase and a low-density phase. By (dynamic) phases in 
this paper we mean not a particular distribution P[a] of arrows a, but a particular 
type of evolution of those distributions Pt[a\. 

To describe the condensed phase let us arrange all the arrows (except the one at 
the root of the tree) downwards, along the edges of the tree. In this case the walker 
starting at the root at time To = returns to the root at time instants T\ = 2, 
T2 = 8, T3 = 22, .... In general, the f th return to the root takes place at the time 
instant T f = 2^+ 2 - 2f - 4. 

There is a growing domain of visited sites — the explored domain — which 
penetrates the k th generation of the tree at the time instant tf. = 2 k+1 — k — 2, 
k — 1, 2, At the time instant = 2 fc+2 — 3k — 3 the domain swallows the A; th 
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generation completely, and the walker heads toward the root. 

If we denote g max (t) the highest generation visited by the walker by the time t, 
then the formula for t k yields 

log 2 (t) - 1 < g max (t) < log 2 (*), for t > 4. 

Analogously, if we denote g c (t) the number of generations completely explored by 
the time t, then the formula for s fc yields 

log 2 (t) - 2 < g c (t) < log 2 (0 - 1, for t > 3. 

Thus, for the downward initial arrangement of the arrows, the growing explored 
dbTQdb IS db "compact" domain of the tree. The height of the domain (measured in 
generations) grows with time as log 2 t. Below the highest visited generation the 
density of visited sites is 1, above that generation the density of visited sites is, of 
course, 0. This is the condensed phase of the Euler walk. 

Another regime - - the low-density phase — is obtained if we begin with the 
upward (left or right) initial orientation of the arrows. In this case the walker goes 
straight toward the top of the tree. The density of visited sites in the f th generation 
at time t is 2~^ +1 (for / < t), which tends to with / justifying the name the 
low-density phase. Of course, once the walker reaches the top of the (finite) tree 
it turns back and gradually stomps the whole graph. Therefore (as it should be) a 
clear-cut distinction between the two phases exists only on an infinite tree. 

For a random initial arrangements of the arrows we obtain a phase which is 
a perturbation of either the condensed or the low-density phase. As we will see 
in the following sections, the transition between the two phases takes place when 
2Pr[\] + Pr[/] = 1. 

The rest of the paper is organized as follows. In Section 2 we investigate the 
properties of the condensed phase: the returns of the walker to the root, and the 
growth of the explored domain when 2Pr[\] + Pr[/*] < 1. In Section 3 we show 
the absence of a compact domain of visited sites if 2Pr[\] + Pr[/"] > I. In Section 
4 we repeat the program of Section 2 at the critical point 2Pr[\] + Pr[/"] = I. 
Traditionally, the last section is devoted to a discussion of the results obtained in 
the previous sections. 

2 The condensed phase. 

Let the internal arrows be initially arranged independently of one another, and 
according to the distribution Pr[\] = p, Pr[/*] = q, Pr[|] = 1 —p — q. One can map 
every initial configuration of arrows into a realization of a discrete-time branching 
process according to the following rules. 

Place a particle at the root of the tree. This particle produces exactly one 
descendant — a particle which is placed at the site of the first generation of the tree. 
From the first generation on, a particle produces either 0, or 1, or 2 descendants 
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Figure 3: An initial arrangement of the arrows, the corresponding first return to 
the root of the Euler walker (solid lines), the first-return cluster of the associated 
branching process (discs), and the buds (spades). At the next visit to a site with 
buds an independent first-return cluster will grow from every bud. 

depending on the initial direction of the arrow at the site occupied by the particle. 
If the arrow points downward, then the particle does not have descendants. If the 
arrow points up and right (like the arrow at the first-generation site on Fig. 3), then 
the particle has exactly one descendant placed at the adjacent site in the direction 
of the arrow. Finally, if the arrow points up and left, then there are exactly two 
descendants placed at the two adjacent sites above, see Fig. 3. 

The relevance of the branching process to our main problem stems from the 
following fact. If the branching process degenerates, then the walker returns to 
the root at a finite time-instant T\ equal twice the number of descendants in the 
branching process (not counting the original particle at the root). The first-return 
path encircles the particles in all generations of the branching process, which we call 
below the first-return cluster. 

The above correspondence between paths of the walker on Cayley tree and re- 
alizations of the branching process allows one to employ the elegant technique of 
generating functions and the main results from the theory of branching processes 
P [5]. First of all recall that if a particle produces k descendants with probabil- 
ity pk, then the branching process degenerates with probability 1 if and only if 
SjfcLi kpk < 1- Hence, the time of the first return is finite with probability 1 if and 
only if q + 2p < 1. The critical case q + 2p = 1 requires a special consideration, 
therefore, in this section we consider only the case q + 2p < 1. 

Lemma 1. Let q + 2p < 1, then the walker returns to the root for the first time 
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at an almost surely finite even time- instant Ti, such that 

mi = ET X - 1 
VarTi = - 



l-(q + 2pY 
4(1 -q) 



l-(q + 2p)f l-(q + 2p) 



Pr[Ti = 2k] „ ^ ,-3 /2 ( 9 + 2v ^r^))' , 

as — >• oo. 

Proof. Denote X the number of descendants for a particle outside the root of 
the tree. The probability generating function of X is given by 

9(y) = Ey x = l-q-p + qy + py 2 . (1) 

Denote Z the total number of descendants in the associated branching process. The 
probability generating function of Z, f(x) = Ex z , is a solution of the equation, see 

f(x) = xg(f(x)). 



Hence 



/(*) 



2px 



1 — qx — y (1 — gs) 2 — 4p(l — p — q)x 2 . 
Differentiating /(x) and taking into account T\ = 2Z, we obtain 

2 „ 4(1 -g) 4 



mi = £^Ti 



1 - (q + 2p) 



, Var T\ 



(l-(g + 2p))3 l-( g + 2p )' 



The above generating function j{x) = Y^k=oPk xk often appears in the literature on 
branching processes, see, e.g., the paper [7] by Otter. In particular, it is shown in 
that paper that the large-/c asymptotics for p^ = P[Z = k] is given by 



Pk 



\ 



gyl 1 -p- q)/p + 2(1- p- q) 
Airp 



k- 3/2 (q + 2Jp{l-p-q) 



(2) 



which is the announced formula for P[Ti = 2k] in the statement of this lemma. 



Remark 1. Denote Z^ the number of particles of the associated branching pro- 
cess in the k th generation of the tree. The random variable X is the number of 
descendants produced by a single particle as in the proof of Lemma 1. Then the 
distribution of the height of the first-return path, Hi, is given by 

Pr[Hi — k] — Pv[Z k+ i = 0] - Pr[Z k = 0]. 
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It is shown in the book by Harris [6] that the large- k asymptotics of Pr[Zk = 0] is 
given by 

Pv[Z k = 0}~l- Cl (EX) k , 

if EX < 1, where c\ is an unknown positive constant. 

Hence in our case the distribution of Hi decays exponentially with k, 

Pr[#! = k] ~ c(q + 2p) k . 

■ 

Thus, during the first stage of exploration of the Cayley tree (0 < t < T\) the 
walker stomps a first-return path with statistical properties described in Lemma 
1. To visualize the motion of the walker after the first return to the root one can 
imagine that, whenever a site is visited for the first time and X descendants are 
produced in the associated branching process, the walker attaches 2 — X buds to 
the site, see Fig. 3. During the second stage of exploration (after the first return 
but before the second return to the root) the walker follows the first-return path, 
but, whenever a bud is encountered, it wonders off the beaten track and appends to 
the existing path a new circuit, which (unless hitting the boundary) is statistically 
equivalent to the first-return path, see Fig. 4. 

Lemma 2. Let the walker return to the root for the first time at time T]_. Then 
the first-return path has exactly 1 + -Ti attached buds. 

Proof. Recall the following standard representation for the number of descen- 
dants, Zfr, in generations k = 2, 3, . . . of the associated branching process 

z, = xf\ 

z, = x? + x? + ... + xtl 
z 4 = x{ 3) + xP + ... + xg, 

and so on, 

where Xj® is the number of descendants produced by the k th particle from the 
/th g enera tion. All the random variables X® are independent and have the same 
distribution as the random variable X. Note also that Z\ = 1, and Z n+ \ = 
whenever Z n = 0. 

Then we have the following formulae for the number of buds bk, in generations 
k = 1,2,3,... 

h = 2 - Xj 1} = 2 - Z 2 , 

b 2 = 2- X[ 2) + 2 - Xf ] + . . . + 2 - Xf 2 = 2Z 2 - Z 3 , 

b 3 = 2 - X{ 3) + 2 - X { 2 ] + . . . + 2 - Xf 3 = 2Z 3 - Z 4 , 
and so on. 

Since for q + 2p < 1 only a finite number of Z k have non-zero values, the total 
number of buds on the first-return path is given by 

oo oo oo 

Bi = E h = Y,( 2Z k-i - z k) = i + E z k- 

k=l k=2 k=l 



7 



1 

Figure 4: The initial arrangement of the arrows, the corresponding second return to 
the root for the Euler walker (solid lines), and a new set of buds (spades). 

The total number of descendants in all generations is ^T±, hence B± — 1 + \T\. 



Theorem 1. Let q + 2p < 1, then the Euler walker returns to the root infinitely 
often at (almost surely finite) time instants Ti,T2,T%, . . .. Moreover, the sequence 
of normalized differences 

Yn= (1 + \ET,Y ' n=1 > 2 > 3 >---> 

is a positive and uniformly integrable martingale, E\Y n \Y n _i, . . . , Y"i] = Y n -\. 

Proof. In order to return to the root for the second time the walker has to repeat 
the first-return path and to create new first-return circuits at each of the B\ buds. 
Hence for the time of the second return to the root we obtain 

(2) 



T 2 -T 1 = T 1 -T + t[ 2) + r 2 (2) + . . . + rg 



(2) 

where To = 0, and rj are independent random variables with the same distribution 
as the first-return time T\. A verbatim repetition of the argument from the proof 
of Lemma 2 shows that on each of the new circuits attached to the first-return path 
the walker creates 1 + \rf^ buds, j = 1,2,..., B\. Hence the total number of buds 
on the second-return path is given by 

.7=1 
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We have essentially the same scenario for any return to the root. For the time 
of the n th return to the root we obtain 



T n T n _i T n _ 



T 



1 - l-n-2 i ^ ■ j 
J'=l 



E 



» 



The number of buds on the n th return path is given by 



B n = + - E 



» 



2 ^ J 



Since _Bi = 1 + ^Ti, we can rewrite the last equation as 



(3) 



Bn= Et'- 



(4) 



where = 1 + |rj n ' ) are independent random variables with the same distribution 
as Bi. 



Induction and the obtained relationships for T n and B n yield 
B n = 1 + \{T n - T„_i), for any n > 1. 



(5) 



Indeed, Lemma 2 says that in the case n = 1 this formula is correct. Suppose that 
the formula is also correct for n = k. Then the relationships for B n and T n yield 



B k+1 = B k + ^Yl ? } 



1 + \{Tk — Tk-i) + \(T k+ i — T k — T k + T k _x) 
1 + \ {T k+ i — T k ). 



Hence B n — 1 + |(^n — ^n-i) is also correct for n — k + 1, which completes the 
induction. 

Now one can calculate the following conditional expectation 



E [T n+1 - T n + 2|T n - T n _x + 2] = T n - T n _x + 2 + E 



E 



_(n+l) 



T n — T n _i + 2 



Hence the sequence 



= (l + ±ET 1 )(T n -T n _ 1 + 2). 



^ri = — T. i ~ \ 5 71 — 1,2,3,..., 



is a positive martingale. 

Since EY n = EY\ = 2, we have Pr[F„ < oo] = 1, for any n, which implies the 
almost sure fmiteness of the return times T n , n — 1, 2, 
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The relationship sup n E{Y%) < oo is a well-known sufficient condition for the 
uniform integrability of the sequence {Y n }^ =1 , see, e.g., the book by Shiryaev [TT] . 
In our case Eq. Q yields 



■Sji+i = E{T n+ i — T n + 2) 2 

= E(T n - T n _i + 2) 2 + 2E 



3=1 



E\Y,rr +1) 



On calculating the expected values with the help of the tower property we obtain 
the following simple recurrent relationship 

s n+ i = s n (l + \ET x f + Var(T 1 )(l + \ET x ) n . 

Solving the recurrent relationship we obtain 



Sn+l 



2 (2 + f|) (1 4- \ET^ - 2^|l> (l + \E Tl )" 



(6) 



Hence 



sup 



1 + lETi 



2n 



< OO, 



implying the uniform integrability of the martingale Y n , n — 1, 2, . . .. 



Corollary 1. Let q + 2p < 1, then for almost all initial arrangements of the arrows 



T n — T n _i + 2 
Jim — 5 r — = Y, 



where Y is a random variable with a proper distribution (Pr[y < 00] = 1). The 
expected value and the variance of the random variable Y are given by 

EY = 2, Var(y) 



ET 1 {2 + ET X )' 

Proof. Since the sequence {Y n }^L 1 is a positive martingale, the Doob martingale 
convergence theorem, see, e.g., the book by Shiryaev [TT], tells us that lim^oo y n = 
Y, where Y is a random variable with a proper distribution. Since the sequence 
{Ynj'^Li is uniformly integrable EY = lim^oo EY n = 2. 
Equation yields 

,. , 4Var(T!) 

lim EYl = 4 + 



ET X (2 + ETi) 



To show that .EF 2 = lim^^ EY^, we need the uniform integrability of the se- 
quence {Y^}^ =1 . To that end one can use the sufficient condition sup n EY^ = 
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sup„ E(Y 2 ) 3 / 2 < oo. One can check by a direct calculation similar to that used in 
the proof of Theorem 1 that the sufficient condition is indeed satisfied. Hence 

Var(F) = EY 2 - (EY) 2 4V&r(Tl) 



ET\{2 + ETi) ' 



Corollary 2. Let q + 2p < 1, then 

2 + ETi r 



ET n 



ETi 



(l + \ETi 
{\ETi) 



- 1 

2n+l 



- 2n, 



lim 



T 



1 + I^Ti 



n+1 



2F 



almost surely, 



where the random variable Y is identical to the one from Corollary 1. 
Proof. Recall that the martingale {Y^}^ is defined by 



Y, 



Tk — 7fc-i + 2 
(1 + \ETi) k ' 



(7) 



Taking the denominator to the l.h.s. and summing over k from 1 ton one obtains 



T n = J2^ + lETi) Y k -2n 



k=l 



where we have used T = 0. Since EYi = 2, a summation of the geometric series 
yields 

2 + ET r / 1 \ n i 
ET n = 2^ f - 1 [(l + \ETi) -l]-2n. 

Equation (jSj), the martingale property £(Yfc|Y/) = Y), for / < fc, and straight- 
forward calculations yield the main asymptotics of the variance Var(T n ), Eq. (|7j). 

Since lim^oo Yk = Y (almost surely), an application of a standard technique 
from analysis to Eq. (jSJ) yields 



lim 

n— >oo 



T 



(l + \ETi 



fe=l 



fc-n-1 2Y 

Yu = 

ETi 



The last Corollary describes in detail the large-n behaviour of the n return 
time T n in the subcritical regime, where ET\ < oo. The following crude bound will 
be helpful at the critical point. 
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Corollary 3. 



T n — T n _! < T n < 2(T n — T n _i). 



(9) 



Proof. For the number of buds on the n th return (to the root) path we have 
B n > 2_B„_i, hence B n -\ < 2~ l B n . Summing Eq. ([51) we obtain 



Theorem 1 and its corollaries give a fairly comprehensive description of the fre- 
quency of return to the root. Our next aim is a description of the height of the 
domain of visited sites. Remark 1 describes the distribution of the highest visited 
generation at time T\. Investigation of the height of the domain at later times is 
a much more delicate problem. We will find the asymptotic behaviour of the den- 
sity of visited sites, Vk(t), in the k th generation of the tree, defined as the ratio of 
the number of sites visited by time t to the total number of sites, 2 k ~ l , in the k th 
generation. 

In order to describe the growth of the domain of visited sites on the Cayley 
tree, let us consider an arbitrary branch w n = (ex, . . . , e n ) of the tree, where ei, 
I = 1,2, ... ,n are the segments (edges) of the branch, see Fig. 5. With any edge ei 
one can associate an "energy" ei as follows. The energy of a left edge ej (like e a on 
Fig. 5) is equal to 0, if the arrow at the bottom of the edge e/ points along the edge, 
and Ek = 1 otherwise. The energy of a right edge e/ (like on Fig. 5) is equal to 
1, if the arrow at the bottom of the edge t\ points down, and E\ = otherwise. In 
other words, the energy of an edge e/ is equal to 1, if the arrow at the bottom of the 
edge causes the walker to deviate from the Euler circuit, and the energy is equal to 
if the walker passes the edge "effortlessly" . The energy of a branch w n is the sum 
of the energies of its edges. 

The domain of visited sites swallows up the edges of a path w n as follows. During 
the time interval [0, T\\ (before the first return to the root) the domain swallows all 
the edges of the path w n till the first obstacle — the first edge ei with E\ = 1. 
During the time interval [T!,T 2 ] (after the first return but before the second return 
to the root) the domain of visited sites swallows up the edge e/ and all zero-energy 
edges which follow ei until the second obstacle — the second edge e m with non-zero 
energy, and so on. During the time interval [Tj, Tj + i] (after the j th return but before 
the j + 1 th return to the root) the domain of visited sites swallows up all the edges 
between the j th and j + 1 th edges with non-zero energy. Thus, the number of visited 
sites in the k th generation at time T m is equal to the number of paths Wk with less 
than m obstacles, or, equivalently, with the path energies E{wk) = S/ :e , ew . Si less 




Therefore T n < 4B n — 2n, and using Eq. <^ again we obtain 



T n — T n _i < T n < 2{T n — T„,-i). 



than m. 
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Figure 5: A branch u? 5 (path) of the Cayley tree, its edges (ej, . . . , es), and the 
associated random "energies" (ei,E2, ■ ■ ■ ,£5)- The energies of edges growing from 
the same site of the tree, like e a and £&, are not independent. 



Let us consider the following sum (partition function) 

efc = $>x P [-/3£K)], 

where the summation runs over all branches w k of a tree with k generations. We 
have 



@fc = J2 #i w k ■ E(w k ) = n}exp[-Pn]. 

n=Q 

Hence, the large k limit of k~ l \n® k is the Legendre-transform of 

u(y) = lim fc -1 ln#{w ft : E(w k ) = [ky]}, 

k— >oo 

where [ky] is the integer part of ky. 

On the other hand, the sum n is almost identical to the partition function of 
a directed polymer on a Cayley tree, see [3j. The difference between 6 n and the 
partition function in [3] is that not all the energies E\ are independent. Indeed if two 
edges e a and e& grow from the same site of the tree, see Fig. 5, then 

Pr[e a = 1, e h = 1] = 1 - p - q, Pr[e a = 0, s b = 1] = 0, 
Pi[£ a = 1, e b = 0] = g, and Pr[e: a = 0, e fe = 0] = p. 

Nevertheless, the large-A; asymptotics of k~ l In 0*. can be found by virtually verbatim 
repetition of the derivation from [3]. In particular, if we denote A k the cr-algebra 
generated by the random energies of the first k generations of the tree, and define 



[{2-2p-q)e~P + 2p + qf 
13 
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then the stochastic sequence {Mk,Ak}kLi is a positive martingale, and EM^ = 1. 
Using the martingale technique from [3] we obtain. 



Proposition 1. If < 2p + q < 1, then 



f(P)= lim ferine* 

fe^oo 
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(2-2p- q)e-P + 2p + q] , if /3 < fl. 
(2 - 2p - g)e- ft + 2p + g| , if /3 > (3 C 



where /3 C is the positive solution of 



In 



(2 - 2p - q)e~ p + 2p + q 



(5{2p + q)et 



2-2p-q+(2p + q)e< 3 ' 



While if 1 < 2p + q < 2, then 



f(/3) = lim k' 1 In 9 fc = In [(2 - 2p - q)e~ p + 2p + q 

k— »oo L 



(10) 



Lemma 3. The logarithmic asymptotics of the number of path with the energy 
, y 6 (0, 1) is given by 



v{y) = lim & 1 ln#{w fe : £(to fc ) = [%]} 

fc^oo 

' , 2-(2p + g) 2p + g 

2/ In + I 1 -y) ln 



y 



i-y 



(12) 



where [x] + = max(x, 0) is the positive part of x. 

Proof. The free energy f(f3), given by Eqs. ffTUl) and ffTTT) . is the Legendre trans- 
form of the logarithmic asymptotics v{y)- Namely 



Therefore 



f(J3)= max Rfy + *%)]■ 
i/e[o,i] 



i/(y) = mm [/3y + /(/?)] . 



Solving the minimization problem we obtain Eq. (TT2"j) . 

■ 

An inspection of the function z/(y) shows that there are around [2 — (2p + q)] k 
branches Wk containing k obstacles for the walker to overcome. At the same time 
there are a few branches with only around [y*k] obstacles, where y* G (0, 1) is a 
solution of the equation 

yln 2-(2 P + g)+(1 _ y)ln 2p+ i ^ 



y 



i-y 



Hence, there exists a growing with time gap, of the width m(l/y* — 1) generations at 
time T m , between the highest visited generation and the highest completely explored 



14 



generation of the Cayley tree. Therefore, neither generation is likely to be a sensible 
measure of the height of the domain of visited sites. 

It is a common practice in situations like that to concentrate ones attention on 
typical branches of the tree. Therefore, we define the height of the domain of visited 
sites as a number (function) H(t) ~ hint, such that the density of visited sites in 
generation xlnt at time t, v x \ n t(t), tends to zero with t if x > h, and v x \ n t(t) — > 1, 
if x < h. We will see shortly that this definition is a sensible one for the problem 
under consideration. Of course, the choice of the asymptotic form H(t) ~ hint is 
specific to Cayley trees, and was actually made after the density of visited sites was 
calculated. 

The logarithmic asymptotics u{y) attains its maximum, In 2, at y = p + \q. 
Hence, the typical branches w k have the energy E{wk) ~ k{p + |g). Thus, the 
domain of visited sites swallows up a typical branch w k of the tree after k(p + |g) 
returns to the root. 

Theorem 2. Let q + 2p < 1, then the height of the domain of visited sites, H(t), 
grows as logarithm of time, 

H(t) 



(p + \q) In (l + i£Ti) ' 



Proof. As follows from Corollary 2, the number of returns to the root by time t for 
the walker is given by 



hit 

m ~ 



ln{\ + \ET 1 



as t — > oo. The asymptotic number of obstacles in a typical branch w k of the Cayley 
tree is given by k(p+ as k — > oo. Hence, the typical penetration after m returns 
to the root is approximately m/(p+ \q) generations, while the typical penetration 
by time t is 

Hit) ~ t r generations. 

(p + \q) In (l + lET ± ) 

■ 

Unfortunately it is difficult to go beyond the logarithmic asymptotics u{y) of the 
number of paths w k with the energy E(wk) = [ky]. Nevertheless, one can guess that 
the number of paths with the energy E{wk) ~ k(p+ \q) + \fku is controlled entirely 
by the quadratic term in the Taylor expansion for u{y) at y = p + |g. If this is 
indeed the case then, in the spirit of the local limit theorem, we obtain 

# \w k : E(w k ) = k(p + \q) + Vku} ~ exp ku(p + \q) + \v' '(p + \q)u 



2 k c 

Tk exp 



u 2 



2(p+\q)(l-p-\q) 



(13) 



The density of visited sites in generation n at time T m is given by 
v n (T m ) = -± =I 'E#{w n :E(w n )=f}. 

Z f<m 
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Approximating the sum by an integral (very much like in the normal approximation 
to the binomial distribution) and taking into account Eq. (|T3l) one obtains 

1 r[m-n(p+q/2)]/^n ( X 2 \ 

Vn{Tm) ~v^L dxex v{-^)i (14) 

where a 2 = (p + \q){l — p — \q). 

We summarize the above discussion by a hypothesis which might well be true. 

Hypothesis 1. The width of the boundary of the domain of visited sites of size 
n generations grows with n as sjn. The drop of the density of visited sites on the 
boundary from 1 to is described by the error function, see Eq. (|14p . 

■ 

Note that at the critical point 2p + q = 1 the variance a 2 in Eq. ([HI) reaches its 
maximal value, |, but remains finite. Therefore the density profile of the domain of 
visited sites does not disintegrate as we approach the critical point. Instead, as 2p+q 
approaches 1, the walker tends to spend more and more time in long (low-density) 
excursions away from the compact domain of visited sites. Those long excursions 
do not create new compact visited domains, somewhat like water poured into sand 
does not create puddles. 



3 The low-density phase. 

Let now q + 2p > 1. In this case the associated branching process degenerates with 
probability x* which is a solution of the equation x = g(x) less than 1, see [6l [5], 
where the function g(x) is given by Eq. (CQ). That is, x* = (1 — q — p)/p. A routine 
application of the Borel-Cantelli lemma shows that in this case, with probability 1, 
the Euler walker visits the root (and any given generation of the tree) only a finite 
number of times. 

Let k be large enough to guarantee that only one copy of the associated branch- 
ing process — the copy which does not degenerate — has survived until the k th 
generation. Then the number of visited sites in the k th generation, (after the 
last visit of the k th generation) does not exceeds the number of particles in a single 
copy of the associated branching process. Namely, < W(q + 2p) k , where W is 
a random variable with a proper distribution (P[W < oo] = 1). Since q + 2p < 2 
unless p — 1, we have Vk/2 h — > as k — >• oo. That is, the model is in the low-density 
phase when q + 2p > 1. 

The bound < W(q + 2p) k is a gross overestimation of the number of visited 
sites. Most likely does not grow faster than something like a constant times In k. 



4 The critical point. 

In this section we consider the critical case q + 2p = 1. Like in the sub critical 
case q + 2p < 1, the associated branching process degenerates with probability 1 
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if q + 2p = 1. However the branching process becomes critical, and its properties 
differ substantially from those in the subcritical regime. As we shall see shortly, 
the first moments of all relevant random variables are infinite if q + 2p = 1. As a 
consequence, extraction of properties of the random variables from their generating 
functions is no longer straightforward. 

Lemma 4. Let q + 2p = 1, then the walker returns to the root for the first time 
at a finite (almost surely) time- instant T%, such that 



Pr[T\ = 2k] 



1 



2^Wp 



k . £ls k 



oo. 



(15) 



Proof. Analogously to the subcritical case, the probability generating function of 
the total number of descendants, Z, is given by 



/(*) 



1 + 



1 



2px 



1 



x 



l-x)[l - (l-4p)x] 



(16) 



Using Eq. <^ we obtain 

Pr[Ti = 2k] = Pi[Z = k] 



2^/Tip 



k ^ . cis k 



oo. 



Remark 2. The large-/c asymptotics of Pr[Ti = 2k] makes it clear that ETi = oo. 
It is still desirable to have a deterministic measure indicating likely values of the 
first-return time T\. For that purpose one can use the quantiles Qi(x) — solutions 
of the equation Pr[Tj < Qx] = x. The asymptotic formula (fl5|) yields the following 
equation for approximate values of Qi(x) 

—L^ v k~^ 2 = i-x. 

2 V^P k >Q 1 /2 

Replacing the sum by an integral and solving the obtained equation for Q\{x) one 
obtains Qi(x) ~ np (i- x ) 2 • ^ or va l ues °f x close to 1, the precision of the found 
approximation for Q\{x) is quite reasonable. For instance, in the case p = 0.1 it 
gives Qi(f) ~ 102, while the exact value is Qi(|) = 98. 



Remark 3. Like in the subcritical case, see Remark 1, asymptotic properties 
of the distribution of the height of the first- return path, Hi, follow from standard 
results of the theory of branching processes. It is shown in the book by Harris [6] 
that the large-A; asymptotics of PrfZ^ = 0] in the case EX = 1 is given by 

1 - p ^ = °i ~ 
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Hence, the distribution of Hi displays a power-law decay, 

Pr[# a = k] = Pr[Z k+1 = 0] - Pr[Z k = 0] ~ 

■ 

In order to investigate the distribution of the return to the root instants T 2 , T3, . . . 
let us first find the probability generating functions G2(x), G 3 (x), . . . for the number 
of buds B 2 , B 3 , . . . on the corresponding paths. Using Eq. (TjJ and the tower property 
one obtains 

G n {x) = Ex B - = E{xf{x)) B — = G n ^{ip{x)), 

where (p(x) = xf(x) is the generating function of B±, and f(x) is given by Eq. (fT6]) . 
It is clear now that G n (x) is the n th iteration of <p(x), that is, 

G n (x) = <p(<p(. ..(p(x) ...)). 
n times 

Hence G n (x) = ip{G n _i(x)) as well. 

Theorem 3. Let q + 2p — 1, then the walker returns to the root infinitely often 
at (almost surely) finite time instants T\, T2, T3, . . .. Moreover 

Pr[T ra - T n _! = 2k] - 2np(1 _ 2 _ n) 1 pl _ 2 -„ A;1+2 -. as fe ^00. (17) 

Proof. The probability Pr[T n — T n „i = 2k] is given by the integral 

where C is a sufficiently small closed contour encircling 0, and G n {x) is the proba- 
bility generating function of B n . To find the large-fc asymptotics of this integral we 
adapt the contour integration from [7j . For that purpose we have to know analytical 
properties of the generating functions G n (x). 
By definition 

00 

G n {z) = Y J ^[Bn = l]z\ 

1=0 

hence the function G n (z) is analytic inside the unit circle {z : \z\ < 1}. Since 
G n -i(l) = 1, and G n (z) = (p(G n -x(z)), the point z = 1 is a branch point of G n (z). 
Since |G n (e J:r )| < 1 for any real x G (0, 27r), the point z — 1 is the only singularity 
of the function G n (z) on the boundary of the unit circle {z : \z\ < 1}. 

From the explicit formula for the function ^p(z) it is clear that the generating 
function G n (z) has only a finite number of points of non-analyticity. Hence, there 
exists a disc A n = {z : \z\ < a n }, with a n > 1, such that z — 1 is the only singularity 
of the functions Gf.(z), k = 1, 2, . . . , n in A n . 



18 



Denote D n the boundary of the disc A n with a radial cut running outwards from 
x = 1. The generating function G n (z) can be written as follows 



GJz) 



1 - aJl - z) 



+ a-z) 



fn{z), 



where f n (z) is analytic and bounded inside D n : \f n (z)\ < b(p) < oo. Indeed, we 
already know that the function G n (z) is analytic inside D n . Since 



fn(z) 



GJz) - 1 + aJl - z) 



(1-z) 



2 -n + l 



it must be analytic inside D n as well. 

To show that f n (z) is bounded inside D n we can use induction. The function 
fi{z) is obviously bounded in any circle with finite radius. Assume now that f n (z) 
is bounded in any circle with finite radius for n = k, then for n = k + 1 we obtain 



G k+1 (z) = v{G k {z)) = 1 



where 



-a-z) 



fk+x(z 

1 

2^ 



1-1 



a k -(l-z) 2 f k {z) 



a, -(I - :)- './,(;) [1 - (1 - W^Az)]-J^ 



Hence f k+ i(z) is bounded in any circle with finite radius as well, completing the 
induction. 



From the above equations we obtain the recurrent relationship a k +i = y/a^p" 1 , 
with the initial condition a x = y/p^ 1 . The solution of this recurrent relationship is 
given by a k = p~ 1+2 . 

On substitution of Eq. ffTS]) in the integral representation for the probability 
Pr[T„ — T n _! = 2k] we obtain 



Pr[T n - T, 



n-l 



2k] 



[l-zf- n ^ | 1 
27rz Jc z k+2 2i\i Jc 



1-z 



fn(z) 



-D k a r . 



k + 1 



+ 



27ri Jd 



(1 - zf~ n+1 f n (z) 



yk+2 



yk+2 



dz. 



dz 



Since the function f n (z) is bounded inside D n , the remaining integral is of the same 
order as 

:i-x) 2 ~ n+1 , „(2- n+i \ 



x 



k+2 



dx = O 



k 



Therefore 



Pr[T n -T, 



n-l 



2k] 



1 



2 n T(l-2- n )p 1 ~ 2 - n k 1 + 2 - n 
19 



as k — ► oo. 



Finally, note that Pr[T n — T n _i < oo] = 1, and according to Eq. (Q we have 
T n < 2(T n — T n _i). Hence, all return to the root instants T n are almost surely finite. 

■ 

Theorem 4. Let q + 2p = 1, then the median of the height of the domain of 
visited sites grows with time as the iterated logarithm 2 log 2 log 2 1. 

Proof. If q + 2p — 1, then the number of obstacles in a typical branch Wk of a 
Cayley tree is ~ k/2. Therefore it takes ~ k/2 returns to the root for the domain 
of visited sites to reach the k th generation of the tree. 

The median m(n) of the duration of n th return loop satisfies 

oo °° 1 1 

Pr[T n - T n _i = 2k] ~ OnrC\ _ 9-^1-2-™ i.1+2"™ = O" 

k>m(n)/2 k>m(n)/2 ^ 1 K 1 z )P h z 

Replacing the sum by an integral and solving the equation for m(n), we obtain 

m{n) ~ 2 2 "c, 

as n — > oo. That is, with probability |, it takes over 2 2 "c time units for the walker 
to complete the n th return path. 
According to Corollary 3 

T n — T n _i < T n < 2(T n — T ri _i). 

Hence the median of T n is between 2 2 "c and 2 2 ™ +1 c once n is sufficiently large. 



The k th generation of the tree is reached with probability \ at a time t ~ c 2 
Solving the equation t = c2 2k/2 for fc, we obtain 



2 fe/2 



k ~ 2 log 2 log 2 1, 



as t — > oo. 



5 Discussion and concluding remarks. 

The analysis of the previous sections can be generalised to the case of a Cayley tree 
with the branching ratio b > 2 at the expense of extra technical efforts. Let the arrow 
directions at every site be numbered counterclockwise 0, 1, 2, . . . , b starting from the 
direction towards the root. Let also Pr[X = k] = k = 0, 1, . . . , b be the initial 
distribution of arrow directions at every site of the tree. Then we can associate a 
realization of a branching process to every initial configuration of arrows as follows. 
From the first generation on, a particle of the associated branching process at a 
particular site of the tree has k G {0, 1, . . . , b} descendants if the arrow at that site 
points in the direction number k. The new particles are placed immediately above 
the parent at the adjacent sites in the directions 1,2, ... ,k. 
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The associated branching process is critical if Y%=i kpk = 1- Already for b = 3 
the explicit formula for the generating function of the total number of particles in 
the branching process, f(x) = Ex z , becomes very cumbersome. For b > 4 we lose 
the luxury of explicit formulae completely. Nevertheless, the results of Lemma 1 
are not difficult to derive for the case of general b. For the condensed phase the 
formulae for differentiation of implicit functions yield 

9 

m x = ETi 
VarTi 



i-eLiW 

4ELo(* - 1) 2 P* 



The tail of the first-return probability, P[Ti = 2k], can be described in terms of a 
positive solution, x*, of the equation 



k=2 



Namely, see [7j, 



Pr[Ti = 2k] 



/(**) k -3/2 [ft** 



\ 2TTf"(x*) \ X 



as k — > oo. 

We see that the properties of the first-return time for b > 2 are qualitatively 
similar to the analogue results in the case of the branching ratio 2. In a similar way 
all the conclusions of the previous sections can be generalized to the case b > 2, and 
the generalization does not produce a novel behaviour. 

Of course the Euler walk on a Cayley tree is only a toy version of Euler walks on 
2D or 3D lattices. Nevertheless we believe/hope that some of the main features of 
the Euler walk described in this paper are also present in finite-dimensional cases. 
In particular, we believe that finite-dimensional walks also have the condensed and 
the low-density phases, and a transition between them. 

Martingales might prove to be also useful for investigation of the finite-dimensional 
walks, but in what way and to what extent is yet to be discovered. Some general 
properties of the growth of the domain of visited sites on 2D lattices might be similar 
to those found in the present paper. In particular, the drop of density from 1 to in 
2D case might still be described by the error function, cf. Eq. (j!4j) . The relationship 
between the size of domain and fluctuations of its boundary might still be the same 
square- root law as in Hypothesis 1. It is possible to state a few more similar hypoth- 
esis, however, the last one already sounds very bold, and it might be dangerous to 
continue any further. In any case, analytical investigation of the growth of domain 
of visited sites for finite-dimensional lattices looks like a very tough problem indeed. 

Monte Carlo simulations for square lattices with equally likely initial directions 
of arrows at every site were conducted in the papers [10J [9] . The simulations show 
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that for 2-D square lattice the radius of the domain of visited sites, R(t), grows with 
time as R(t) ~ c t 1 / 3 . It was also conjectured that on a 3-D square lattice and in 
higher dimensions we have a diffusive behaviour, R(t) ~ ct 1 / 2 , because the walker 
does not return to the cluster of visited sites frequently enough. The last conjecture 
seems to imply that as the lattice dimension tends to infinity the behaviour of the 
walker does not become more and more similar to that of a walker on a Cay ley tree. 
However, on the basis of results obtained in this paper one can put forward the 
following alternative interpretation of the "diffusive" behaiviour of the walker on a 
3-D lattice. It might be the case that the Monte Carlo simulations for 3-D lattices 
were simply conducted in the low-density phase where a compact domain of visited 
sites is not formed. Changing the initial distribution of arrows one can get into the 
condensed phase, where the radius of the domain of visited sites grows, presumably, 
as R(t) ~ c d t^ d+1 \ converging to the (Cayley tree) logarithmic behaviour as the 
lattice dimension d — > oo. 

Something similar actually happens on Cayley trees as well. If we take a Cayley 
tree with the branching ratio 6 = 2, then the equally likely initial distribution of 
arrows Po = Pi = P2 = I puts the walker at the critical point p\ + 2p 2 = 1. 
While if we increase the branching ratio to 3, then the equally likely distribution 
Po — Pi — P2 = P3 = \ corresponds to the low-density phase p\ + 2p 2 + 3p 4 > 1. 

It was already known that branching processes are relevant to and, in fact, 
provide a mean-field description for some model of self-organized criticality, see, 
e.g. [H[12]. Although branching process are also relevant to Euler walks, the latter 
apparently belong to a somewhat different class of models, since instead of fixed 
values for the standard set of critical exponents, we have a whole spectrum of those. 
Indeed, instead of the mean-field exponent r = 3/2, describing the distribution of 
the size of avalanches, we have the sequence r n = 1 + 2~ n , n = 1,2,..., beginning 
with 3/2. 

Due to the infinite memory of the Euler walk it is difficult to calculate the 
moments of the walker's location, Ex k (t). It is a pity, since the second moment of 
the walker's location for the simple random walk on, say, 2-D lattices, Ex 2 (t) = ct, 
is one of the main characteristics of that random process. To partially fill this gap 
we will extract some information on the behaviour of the second moment from the 
results obtained in the previous sections. This information might provide clues for 
explanation of a bizarre behaviour of Ex 2 (t) for certain versions of Euler walk on 
2-D lattices [8]. It is instructive to compare at the same time the behaviour of the 
Euler walk on a Cayley tree and the simple random walk on a 2-D square lattice. 

Both the Euler walk on a Cayley tree in the condensed phase and the 2-D simple 
random walk are recurrent. Here, however, similarities end. While the expected 
return-time (and even the variance) for the Euler walk is finite, the expected return- 
time for 2-D random walk is infinite. As a consequence we have monotonically 
increasing variance of the walker's location for the 2-D random walk, Ex 2 (t) = ct. 
On a Cayley tree the walker returns to the root at time instants T n with ET n < 00, 
n = 1, 2, . . ., see the explicit formulae in Corollary 2. If Var T\ <C [ET\) 2 , then the 
returns to the root in the logarithmic scale take place almost periodically, lnT„ ~ 
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nln(l + ^ETi), as n — > oo. On the other hand, if VarTx ^> (ETi) 2 , then the 
periodicity in the logarithmic scale turns into chaotic behaviour without any visible 
pattern. 

While in the latter case one can not rule out the monotonic increase of Ex 2 (t), 
in the former case one certainly has a nearly periodic vanishing of Ex 2 (e t ). If the 
magnitudes of VarTi and (ETi) 2 are comparable one should have an intermediate 
situation with visible deviations in the shape of Ex 2 (t) from a classic ct 7 behaviour. 
As we approach the critical point q + 2p = 1, the variance 

iter,- 4(1 - q) 



[l-(« + 2j>)]3- 
grows faster than 

(et,) 2 1 



[l-(? + 2p)] 2 ' 

and we lose completely traces of the log-periodic behaviour. 
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